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The field dependence of tfie vortex core size €,{B) is incorporated in tlie London model, in order to 
describe reversible magnetization M{B,T) for a number of materials with large Ginzburg-Landau 
parameter k. The dependence ^{B) is directly related to deviations in M(ln_B) from linear behavior 
prescribed by the standard London model. A simple method to extract (,{B) from the magnetization 
data is proposed. For most materials examined, ^{B) so obtained decreases with increasing field 
and is in qualitative agreement both with behavior extracted from /iSR and small angle neutron 
scattering data and with that predicted theoretically. 

PACS numbers: PACS: 74.50.+r, 74.78. Fk, 74.78.-w 



I. INTRODUCTION 



On the London model 



Despite its simplicity, the London approach is a pow- 
erful tool in describing magnetic properties of the mixed 
state. In fact, short of the full-blown microscopic theory, 
it is the only method available for low temperatures. The 
approach is based on the London equation: 



rn) , 



(1) 



where h{r) is the magnetic field, A is the penetration 
depth (a temperature dependent constant in uniform 
samples) , 4>o is the flux quantum, and r„ are vortex posi- 
tions. For simplicity, the equation is written for isotropic 
materials. This approach fails at distances of the or- 
der of the coherence length ^; still, in materials with 
K — A/^ 3> 1, there is a broad domain of intermedi- 
ate fields (po/X^ ^ i7 <C 4>o/^'^ where the complexity 
of the vortex core contributions to the total energy can 
be disregarded and the London approach suffices for the 
description of macroscopic magnetic properties. 

As far as the equilibrium properties of the flux-line 
lattice are concerned, the pivotal point is the expression 
for the free energy 



F^F = 



In 



erjH, 



c2 



(2) 
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The right-hand side here is the interaction energy of vor- 
tices forming a periodic lattice, B is the magnetic in- 
duction. This expression is obtained by transforming 



the sum of pair-wise interactions of vortices to a sum 
over the reciprocal lattice, see e.g., ReL 1. The sum (or 
the integral over the reciprocal plane k) is logarithmi- 
cally divergent so that a cutoff at A: ~ 1/^ '--^ 1/pc is 
needed {pc is the size of the vortex core). This yields 
ln((/)o/27r^2i?). The parameter rj of the order unity is 
commonly introduced to account for uncertainty of the 
cutoff (along with the uncertainty in the lower limit of 
the integral of the order of inverse intervortex spacing 
^ \/B/4>q); e = 2.718... appears in Eq.(2) for conve- 
nience of not having it in the expression for the magneti- 
zation. Again, the energy in the form of Eq.(2) holds in 
intermediate fields Hd <C <C i?c2, the domain existing 
only in materials with large Ginzburg-Landau parameter 
K (-ffci a-nd Hc2 are the lower and upper critical fields). 
Hence, although the length ^ (or the core size pc) does not 
appear in the London equation (1), it enters the energy 
expression (2) through the cutoff and therefore affects, 
presumably weakly, macroscopic quantities such as the 
magnetization and other properties of the mixed state. 

Significant effort has recently been put in studies of 
the vortex core size Pc, see the review by Sonier and 
references therein. ^ Notably, whatever definition of pc is 
adopted, the low temperature Pc (extracted from the /zSR 
data) decreases with increasing applied magnetic field in 
a number of materials such as NbSe2, VaSi, LuNi2B2C, 
YBa2Cu307_5, and CeRu2; their physical characteris- 
tics have little to do with one other, except that all of 
them have a large GL parameter k = X/^ and exhibit 
large regions of reversible magnetic behavior. One can 
add to this list a heavy fermion compound CeCoIns, for 
which the interpretation of small angle neutron scattering 
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data (SANS) requires a similar behavior of the coherence 
length.^ The dependencies Pc{B) for aU tested materials 
are qualitatively similar: when the field increases toward 
Hc2, Pc{B) decreases roughly as In other words, 

in large fields Pc is roughly proportional to the intervor- 
tex spacing. 

A few qualitative reasons for the core shrinking with 
increasing field have been discussed in literature; see the 
review 2. Perhaps the simplest defines the core bound- 
ary as a position where the divergent London current 
C(?io/87r^A^r of an isolated vortex reaches the depairing 
value, i.e., as r ~ ^. In the mixed state, neighboring 
vortices suppress the circulating current by contributing 
currents of the opposite direction. Hence, the depairing 
value is reached at a shorter distance from the vortex cen- 
ter, and consequently pc should decrease with increasing 
field. 

The vortex core size, pc, is of the order of the coher- 
ence length ^ and, in fact, it is often identified with ^(T). 
Strictly speaking, the latter is defined only at the up- 
per critical field: Hc2 {T) = (/>o/27r^^. Nevertheless the 
length ^ is used to describe the mixed state at fields not 
necessarily close to Hc2- The question of possible field 
dependence of ^ has been considered by one of us for 
the isotropic case.^ It was shown that in the dirty limit 
one can use ^ = ((/)o/27r_ffc2)^^^ at any field within the 
mixed phase; the same is true near the critical tempera- 
ture Tc for any scattering strength. However, in general, 
when the field is reduced below Hc2, the value of £,{B) 
increases, an effect that is profound in clean materials at 
low temperatures. Calculations of Ref. 5 are in accord 
with the /uSR results cited above. In the following, we 
denote as ^c2 the value of £,{B) at B = Hc2 to stress that 
in general ^{B) 7^ for B < Hc2- 

Another common theoretical definition of Pc is based 
on the slope of the order parameter A(r) at the vor- 
tex axis r = 0, normalized either to its value A(oo) 
far from the single vortex or to the value of A(a/2) 
half-way to the nearest neighbor in the mixed state: 
1/pc = A'(0)/A(a/2). Recent microscopic calculations 
of this quantity by Miranovic et al. showed a variety of 
field dependent behaviors of Pc at low temperatures de- 
pending on the scattering strength.^ In particular, this 
work suggests that Pc{B) may have a minimum which, 
however, has not been seen in /xSR experiments. There 
are many different ways to define Pc customized for dif- 
ferent experimental or theoretical needs (see, e.g., the 
discussion in Ref. 5). For the purpose of this paper, these 
differences are irrelevant and we use the terms pc{T,B) 
and ^(T, B) as the same. 

In the following we provide experimental data for the 
field dependence of the reversible magnetization for a 
single crystal YNi2B2C in a broad temperature region 
to demonstrate the known fact: the data cannot be de- 
scribed by the standard London model. We then derive a 
closed form expression for the cutoff ^{B) needed to rep- 
resent correctly the data M (B) with the help of the Lon- 
don model. We show that £,{B) so chosen is qualitatively 



consistent with the field dependence of ^ recorded by /iSR 
and discussed theoretically. Finally, we demonstrate that 
the new model generates a consistent description of the 
magnetization data for a number of unrelated materials 
with large k. The goal of this paper is to demonstrate 
that ^ (B) can be in principle extracted from the magne- 
tization data, a less demanding experimental procedure 
as compared to pSR or SANS. Within our approach, the 
London penetration depth is field independent, whereas 
the field dependence of £, alone suffices to explain the 
data. 



II. EXPERIMENTAL ASPECTS 

A. Sample preparation 

Single crystals of YNi2B2C, LuNi2B2C, and 
Lu(Nii_a;Cox)2B2C were grown out of Ni2B or 
(Nii_j,Cox)2B flux in a manner similar to the growth of 
other borocarbide crystals. As discussed in Ref. 7, Co 
doping serves as a convenient tool to move from clean to 
dirty limit. 

NbSe2 crystals were grown via iodine vapor transport 
technique (REF) and had Tc w 7.1 K and the residual 
resistivity ratio RRR w 40. 

MgB2 single crystals of submillimeter sizes were grown 
by a high pressure technique as described in Ref. 8. 

A single crystal of Y^Si was grown using a Bridgeman 
method, in which a floating zone was created by rf in- 
duction heating. The samples for investigation were cut 
by a wire saw and oriented using Laue method. Samples 
had typical dimension of 2 x 2 x 4 mm. 

B. Magnetic measurements 

The magnetization measurements were performed by 
using several Quantum Design MPMS systems. In a typi- 
cal experiment, a full M [H) loop was recorded and only 
its reversible part, above the irreversibility field, i?irr, 
was used for the analysis, ffj^r was determined as a field 
where ascending and descending branches coincided or 
were sufficiently close (a weak hysteresis). 

The procedure is demonstrated in Fig. 1. The main 
frame shows raw data with clear range of reversible be- 
havior above Hirr indicated by an arrow. The left in- 
set shows the expanded portion of the raw data in the 
vicinity of Hc2- A small paramagnetic background (from 
the sample, the sample holder and perhaps residual flux 
on the crystal surface) is clearly seen as a linear-in-iJ 
contribution. After this contribution is subtracted, we 
obtain the superconducting diamagnetic signal shown in 
the right inset. The Hc2 is indicated by an arrow. Su- 
perconducting transition temperature was measured in a 
small {H = 10 Oe) applied field. In large fields of our in- 
terest, demagnetization effects are weak; in the following 
we do not distinguish between the applied field H and 
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FIG. 1: (color online) An example of M{H) for YNi2B2C 
at T = 9K. The main plot shows both up- and down-field 
scans and the irreversibility field. The upper inset illustrtes 
how the normal state paramagnetic contribution is subtracted 
with the result shown in the lower inset. The upper critical 
fields Hc2 is indicated by an arrow. 
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FIG. 2: (color online) The magnetization M(ln B) for 
YNi2B2C at T = 2, 5, 9, and f2K. The upper critical fields 
Hc2 are the positions of kinks in M(ln B) not shown in the 
figure. The solid curves are obtained by fitting the data to 
Eq. (4) with the fitting parameters Mo (shown in the inset), 
r;, and a. 



the induction B. As explained below, we do not need 
the sample volume in our analysis so that we can use 
on equal footing field dependencies of the magnetic mo- 
ment (in emu) or of the magnetization (in G) ; we use the 
notation M for both quantities. 



III. MODIFIED LONDON MODEL 

The standard London energy (2) gives an equilibrium 
magnetization that is linear in In B in intermediate fields: 



M = 



dF 



dB 327r2A2 



'0 '7 



B 



(3) 



Hence, the standard London model requires a plot of M 
versus In i? to be a straight line. 



A. YNi2B2C 

Fig. 2 shows reversible magnetization for a single crys- 
tal YNi2B2C in fields parallel to the c axis at 2,5,9, and 
12 K. Clearly, the deviations from the London prediction 
increase with decreasing T; at low temperatures M (In B) 
is far from being linear. 

We note that in many materials with k ^ 1, it is 
difficult to distinguish between a narrow Abrikosov do- 
main near with M oc (i?c2 ~ B) and a broad London 
domain where the magnetization depends on the field 
in a slow, nearly logarithmic manner. For this reason 



the Abrikosov part of M{B) is sometimes discarded al- 
together; this amounts to setting = 1 in Eq. (3).^ Of 
course, this cannot be done for materials with k ^ 1. We 
will follow this simplification in our analysis and indicate 
the cases when this cannot be done. 

To formally account for deviations of the data from 
the behavior prescribed by the standard London formula, 
we add to expression (3) two additional terms: const/_B 
(to correct for the low field behavior) and const -B (to 
account for the high field curvature). These are, per- 
haps, the simplest possible modifications one can think 
about. A restriction upon the constants is provided by 
a requirement that M{Hc2) = 0. Hence 



M = -M, 



In 



c2 



, a in 77 + a 

B B ^ ' 'Ha 



(4) 



where Mq = 0o/327r2A2. Certainly, the form (4) is not 
the only possiblity for representing the available data. 
Other forms were suggested in the literature^ ^"^'^ based 
on different theories and assumptions. We stress that 
the expression (4) is just an empirical formula to rep- 
resent the data. We choose it because of its simplicity, 
and because - as is demonstrated below - it is sufficiently 
flexible to represent the magnetization data in a host of 
materials with very different physical properties. 

The solid lines in Fig. 2 are the data fits to Eq. (4). 
The value of the upper critical field for each T is read 
directly from the raw data as explained in Fig. 1. We are 
left with three fit parameters Mg, ry, and a. Two of these 
are shown in the table of Fig. 2. The inset in Fig. 3 shows 
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that the T-dependence of Mq oc 1/A^ is quahtatively con- 
sistent with the behavior of the superfluid density oc A~^. 
The quahty of the fits is good; hence, the empirical form 
(4) can be used to represent the data reasonably accu- 
rately. 

Examining possible modifications of the London model 
to account for the deviations of M(ln B) of Fig. (2) from 
linear behavior, one should bear in mind the difference 
between the roles of two fundamental lengths, A and ^, 
within the London theory. The length A enters Eq. (1) 
which is the basis of the whole approach. On the other 
hand, the length ^ is absent in the London equation 
and enters the energy expression as an uncertain cutoff 
used to mend the inherent shortcoming of the London 
model. Therefore, considering possible modifications of 
this model, one still has some freedom - however limited 
- in working with ^, unlike the case of A. 

Comparing the data of Fig (2) with predictions of the 
standard London model one wonders why the model 
which describes correctly the field of vortices away of 
their cores, fails badly at low temperatures. From the 
point of view of a consistent London description, the only 
suspicious point in deriving the free energy (2) and the 
corresponding magnetization (3) is the cutoff employed 
which generates the term hi{Hc2 / B) . 




FIG. 3: (color online) ^{B) extracted from the data of Fig. 2 
for YNi2B2C with the help of Eq. (8). The upper curve is 
for T = 2K, the lowest is for 12 K. The inset shows Mo{T) 
extracted from the fits of Fig. 2. 



B. London model modified to accommodate {(B) 

Hence, we write the free energy in the form: 

^^32^^"—^' ^=2<^' 

where ^{B) is the field dependent cutoff (the core size). 
Clearly, ^(i?c2) is the standard coherence length associ- 
ated with Hc2, so that H{Hc2) — Hc2- Then, evaluating 
M = -dP/dB, we obtain: 



M = 



i)o l^^vH , B dH 



327r2A2 \ B H dB 



(6) 



The idea of the following manipulation is to find a field 
H that generates the form (4), or in other words, that 
represents the experimental data. After equating (6) and 
(4), one can solve a linear differential equation for H{B) 
with the boundary condition H{Hc2) — Hc2- The result 
in terms oi h = H / Hc2 and h = B /Hc2 reads: 

ln/.= ^ln5-K^i^^^^±^Ki^. (7) 

This corresponds to the normalized cutoff distance (the 
core radius) 

r(S) = ^=5^"/-exp(i^^^^±4(^. (8) 
Cc2 46 



It is readily shown that the slope of ^*(6) at Hc2 is de- 
termined by the parameter 77: 



dc 



.1 = ^"^- 



(9) 



Hence, when the field decreases from iJc2, ^{B) decreases 
for 77 > 1 and increases for 77 < 1. 

Using Eq. (8), we can calculate the normalized cutoff 
^* (6) responsible for deviations of M {B, T) from the stan- 
dard London behavior for YNi2B2C shown in Fig. 2 since 
we have rj and a representing these data sets (note that 
Mo does not enter Eq. (8)). The curves i*{b) for T = 2, 
5, 9, and 12 K calculated with Eq. (8) are shown in Fig. 3. 

It is worth observing that £,{B) so obtained is qualita- 
tively similar to the B dependence of the core size seen in 
/iSR experiments.^ Moreover, it is argued in Ref. 5 that 
the field dependence of ^ should weaken with increasing 
temperature in accord with Fig. 3. 

Another point to make is only a moderate variation of 
^ which is needed to account for strong deviations from 
the linear M(lni?). For example, at 2K, ^ changes only 
by a factor of 2 over most of the mixed state field domain. 
Since the cutoff enters the energy (5) under the log-sign, 
it might be surprising that such a difference suffices to 
cause a drastic deviation of the 2 K curve in Fig. 2 from a 
straight line. The puzzle is resolved if one observes that 
the field dependence of ^ translates to non-logarithmic 
corrections to the standard London magnetization, see 
Eq. (6). 

The same analysis has been applied to the magne- 
tization data for a crystal of LuNi2B2C [whose crystal 
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FIG. 4: (color online) (a) M{E) for LuNi2B2C for T=7,10 
and 12 K. (b) C(-B) corresponding to the graphs of the upper 
panel. 




FIG. 5: (color online) (a) M(B) for Lu(Nii_^Co,;)2B2C with 
X = 0, 3, and 6%. (b) £,(B) corresponding to graphs of the 
upper panel. 



structure and superconductive properties are similar to 
YNi2B2C], yielding similar results as shown in Fig. 4. 



C. Lu(Nil-,COa;)2B2C 

As mentioned, Ref. 5 argues that the field depen- 
dence of the core size is weakened by increasing tem- 
perature and scattering. In order to study the scatter- 
ing dependence of ^(-B) we turn to a series of crystals 
Lu(Nii_a;Coa;)2B2C in which the mean-free path is pro- 
gressively reduced by increasing the Co content.^ 

In the upper panel of Fig. 5, M{B) is shown for x — 0, 
3, and 6% (for which Hc2 has been independently mea- 
sured) at the same temperature of 2 K. The fit parameters 
T] and a are also shown and the calculated field dependent 
core sizes are given in the lower panel. 

We note that for 6% Co, the ratio of the zero-T co- 



herence length to the mean-free path has been estimated 
in Ref. 7 as exceeding 10, which places this sample close 
to the dirty limit. The core size for this crystal is seen 
to vary only by about 10%, which is in accord with the 
theoretical finding that field dependence of ^ disappears 
in the dirty limit. '^'^ 



D. NbSea 

The superconducting anisotropy of this material is 
stronger than in borocarbides discussed above. As is seen 
in Fig. 6, deviations of M{lnB) from the standard Lon- 
don linearity are profound along with the corresponding 
field dependence of 
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FIG. 6: (color online) (a) M(B) for NbSea. (b) ^(B) corre- 
sponding to graphs of the upper panel. 



FIG. 7: (color online) (a) MgB2, T = 4.6,15 and 25K.(b) 
corresponding S,{B). 



E. MgB2 

The upper panel of Fig. 7 shows A/(lni?) in fields par- 
allel to the c axis of a single crystal MgB2. One readily 
sees a qualitative difference from the preceding examples: 
the curvature of A/(lni?) for T ~ 4.6 K being negative in 
low fields becomes positive in large fields. Still, we can 
fit well the data for all T's to the form of Eq. (4) with 
parameters rj and a given in the table. 

The most interesting feature is that the obtained values 
of rj exceed unity. According to Eq. (9) this means that 
starting from Hc2, C decreases with decreasing field. This 
is shown in the lower panel of Fig. 7. We attribute this 
unusual behavior to the two-gap nature of this material: 
the small gap on the 7r-sheet of the Fermi surface opens 
up in decreasing fields thus causing a decrease of 

If indeed the unusual behavior of ^{B) for MgB2 is 
due to suppression of the small gap in fields of few kG 
for H\\c, and if the suppression field is isotropic^ then 



by going to other field orientation away of the c axis, 
we can push the effect of the small gap out of the high 
field domain of our interest. To check this hypothesis we 
acquired the data for the applied field at 45° to the c axis 
where -ffc2(0) is accessible with our equipment. Figure 
8 shows the result similar to that for Y and Lu-based 
borocarbides. This suggests that, e.g., for T — 10 K, in 
the field domain examined (from Hc2 ~ 4.15T down to 
about 0.4 T or b w 0.1) the small gap is not yet fully 
formed. 

Of course, this interpretation is much too simple 
because for other than -ff 1 1 c orientation the strong 
anisotropy of ^ should be taken into account, the sub- 
ject of our future work It should also be noted here that 
the macroscopic phenomenology of magnetic properties 
of MgB2 is still debated. Despite the two-band nature 
of this material, within the London approach, we employ 
one penetration depth A and one cutoff length ^ when 
describing the vortex lattice in reciprocal space. Judging 
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FIG. 8: (color online) (a) Magnetization of the MgB2 crystal 
in applied field at 45° to the c— axis, (b) corresponding co- 
herence length. Note the difference in the behavior of ^*{b) 
from the case of the field along c of Fig. 7. 



by literature, this point of view is not universally shared 
by all MgB2 community. 



F. VsSi 

Given the slope dHc2/dT = 19.4T/K and = 16.6 K, 
it is likely that the low temperature upper critical field 
of this material exceeds 20 T; no direct measurements of 
Hc2 were conducted in this field range. One might treat 
Hc2 as an extra fitting parameter to be extracted from 
the data on M{B). However, the numerical procedure of 
extracting both Hc2 and rj from the magnetization data 
is unstable because their product enters the formulas we 
use. For this reason we consider here only M{B) for 
T > 13K for which Hc2 was measured. 

It is worth noting that Eq. (4) is good enough even for 
a quite unusual shape of AI{\nB) in this material: the 
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FIG. 9: (color online) Magnetization of the VaSi single crystal 
at r= 13,14, 15 and 16 K. 



curvature of Af(lnB) changes sign in all data we have 
examined. With the help of Eq. (4) we readily find that 
the inflection point is at hi = \J ajia ^-\Y^r\). As is seen 
from the table of Fig. 9, for all curves given, a <C Iut; 
which leads to 



In hi 



a 



' ln?7 



(10) 



Hence the inflection points of M(ln6) are approx- 
imately in the same place, Inbi « 0.5, for all 
temperatures. "'^^ Therefore, the curves M(lni3) are 
concave for b < hi k, 0.16, i.e., at fields hardly in the 
domain of applicability of our high-field model. On 
the other side of the inflection point where the curves 
M{h\B) are convex, we may deal with the Abrikosov 
domain where M is linear in Hc2 — B; it is easy to see 
then that M(lni3) should be convex. Application of 
the London approach in this domain cannot be justified. 
Therefore, there is no point in trying to extract ^(B) 
from the data on VaSi (formally, since all ry's in the table 
of Fig. 9 exceed unity, this extraction would have given 
£,{B) decreasing with decreasing field which would have 
contradicted the /iSR data of Ref. 18). 



IV. DISCUSSION 

The main point of this work is to argue that the 
field dependence of the core size is a generic low- 
temperature property of all sufficiently clean supercon- 
ductors. Moreover, incorporating this dependence in the 
London approach broadens considerably its applicabil- 
ity for describing macroscopic reversible magnetic prop- 
erties. Still, the empirical approach adopted here lacks 
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microscopic justification. If proven correct, the modified 
London scheme calls for revisiting many results obtained 
within the standard London model in which the field in- 
dependent core size or the cutoff are involved. Field de- 
pendencies of the flux-flow resistivity, of the specific heat 
in the mixed state, or the core pinning arc some exam- 
ples. 

One can foresee a number of difficulties in trying to 
develop such a justification. The cutoff size we extract 
from M{B) data is not necessarily the same as the core 
size defined as being proportional to the slope of the or- 
der parameter at the vortex axis: approaching the core 
from the outside to determine the cutoff we may have 
a different result than when examining the core struc- 
ture starting from the core center. Having this in mind, 
it is not surprising that the microscopic calculations of 
Pc oc {dA/dr\r^o)^^ in Ref. 6 do not agree in detail 
with our empirical results (to our knowledge, this is the 
only calculation of this sort that has been tried to date). 
Fig. 10 shows the output of this calculation normalized in 
the manner of this paper for a model isotropic material at 
T — 0.1 Tc and for a few values of the scattering param- 
eter ^o/£ (^0 = hvp/TT^iO) is the zero-T BCS coherence 
length and £ is the mean-free path for the non-magnetic 
scattering) . We note that while the curves generated are 
mostly qualitatively similar to what we extract from the 
magnetization data, this calculation does not confirm our 
assertion about weakening of the field dependence of the 
core size with increasing scattering. 
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cally uniform within the core for k ^ 1, one can use the 
Helfand-Werthamer^^ linearization technique employed 
for calculation of Hc2 also for the core size in the high- 
field mixed state. Within this approach, the order pa- 
rameter near the vortex center satisfies a linear equation 
-^2n2A = A, where 11 = V + 27rA/0o and ^ is found 
solving the basic BCS self-consistency equation. This 
produces ^{T,£;B) in qualitative agreement with what 
we extract from the magnetization data in this work 
(for all cases other than V^Si and MgB2 in field along 
the c axis); in particular, the analytical field dependence 
of ^ obtained in this way disappears when T Tc or 

In our view, the question still remains which of these 
two theoretical approaches, Ref. 6 or Ref. 5, describes bet- 
ter various data on S.{B). An important role in resolving 
the question belongs to experimental studies of how the 
field dependence of the core size affects other physical 
properties of the mixed state. 



A. Flux-flow resistivity 

As an example of such a cross-examination we took 
data on the fiux-fiow resistivity pf from measurements 
of the microwave surface impedance of YNi2B2C.^" The 
data show large deviations of the measured p / from the 
Bardeen-Stephen linear field dependence, pf = pnB/Hc2- 
This formula is obtained assuming a field independent 
core size ^ = ^^2 = ^0o/27riJc2- Clearly, if ^ does depend 
on the field, one has: 
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FIG. 10: (color online) The field dependence of ^ calculated 
microscopically for T /Tc = 0.1 and for a few scattering pa- 
rameters ^0 /£ where is zero-T BCS coherence length and £ 
is the mean-free path for scattering on non-magnetic impuri- 
ties. 



A different approach to evaluation of the core size is 
chosen in Ref. 5: it is argued on physical grounds that 
since A — s- at the vortex center and the field is practi- 
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Hence, for each data set pf{B), we can extract 



m _ pf{B) 



c2 



Pn b 



(11) 



(12) 



in other words, we can delegate deviations from pf oc B 
to the field dependence of ^ and see whether or not the 
obtained ^{B) agrees with that extracted from magneti- 
zation data. 

Utilising the flux-flow resistivity data of Fig. 3 of 
Ref. 20 and applying Eq. (12), we obtain the result shown 
in the panel (a) of Fig. 11. Comparing it with our Fig. 3 
for the same material, we obtain reasonable agreement, 
notwithstanding the usage of different samples in these 
two experiments. 

As another example, two selections of data for 
Y(Nii_:rPt2;)2B2C from Fig. 3 of Ref. 21 were used to 
calculate ^{B)/£^c2, as shown Fig. 11(b). We see that 
in addition to the expected decrease with field, £,{B) is 
suppressed by increasing scattering (i.e., increasing im- 
purity content of Ft) again in accord with the examples 
discussed above. 

Panel (c) of Fig. 11 shows the result of the same ex- 
ercise with a d-wave material (an overdoped crystal of 
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Bi-2201).^^ This example supports the idea that the field 
dependence of the London cutoff is a generic feature of 
type-II superconductors with no direct relation to the 
order parameter symmetry. 

The last panel of Fig. 11 presents the cutoff ^(5) /^c2 for 
two field orientations of MgB2 extracted from the flux- 
flow resistivity data of Ref. 23, again using Eq. (12). For 
the field along ab, both M{B) and Pf{B) data yield qual- 
itatively similar results. This, however, is not the case 
for the field along c as evident by comparing Figs. 7 and 
11(d). Thus it appears that for the two-gap MgB2 with 
its greater complexity, the simple scheme of incorporat- 
ing a field-dependent cutoff to the London model does 
not work for all field orientations. 
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FIG. 11: (color online) ^{b)/^c2 versus b = B/Hc2 for 
Y(Nii-:cPtx)2B2C with x ^ and 0.2 extracted from the 
flux-flow resistivity data of Ref. 21 as explained in the text. 



B. On nonlocality 

Deviations in M(ln B) from the standard London lin- 
ear behavior have been thought to come from the effects 
of the nonlocal relation between current and the vector 
potential inherent for superconductors.^^ The nonlocal 
corrections to London theory turned out to be an effec- 
tive tool in describing evolution and transitions in vortex 
lattice structures. However, it is difficult at this stage 
to sort out what part of the deviations of M from linear- 
in-ln B behavior arises solely from the core-size field de- 
pendence and what part should be relegated to the nonlo- 
cality. In particular, the difficulty comes from analogous 
weakening of the two effects with increasing tempera- 
ture and scattering strength. Looking back to the overall 
satisfactory data analysis of M{lnB) for Lu(Nii_a;Coa;)2 
B2C of Ref. 25 based on nonlocal effects, we note that the 
analysis produced an excessively rapid reduction of the 



"nonlocality range" with temperature for samples with 
elevated impurity content. We have to conclude that 
more high precision experimental results are needed be- 
fore a conclusive judgement is made on the relative im- 
portance of contributions to M{B) due to the nonlocality 
and to the field dependence of the core size. 



C. Why \(B) should be used with caution within 
London theory 

As mentioned above, the basic London equation (1) 
along with the energy expression (2) for intermediate 
fields (or Eq. (5) with an unspecified cutoff) imply that 
the London penetration depth A is a temperature depen- 
dent constant in a homogeneous material. Assuming a 
field dependent A would have changed the London equa- 
tion (1) per se: the quantity A^(/i) cannot be taken out of 
differentiation operators. As a result, the Fourier com- 
ponents of the solution h{k) for a single vortex would 
have been different from 4>o/{l + X^k^) and the energy 
(2) would have been different as well. Therefore, unlike 
the case of the cutoff ^, relaxing the requirement of a 
constant A causes basic changes in the London approach 
and therefore should not be taken lightly. It also worth 
recalling that a constant A is derived from the micro- 
scopic BSC theory (as fc — > limit of the BCS kernel in 
the nonlocal connection between the persistent current 
and the vector potential). To our knowledge, there is no 
microscopic justification for a field-dependent A (in non- 
magnetic superconductors). In other words, the London 
theory is rigid with respect to a constancy of A, unlike 
the case of ^. 

Yet, quite often, analyzing data with the help of the 
London model (i.e., starting with a constant A) it is 
concluded that A should be field dependent. Numerous 
examples are found in the literature on /iSR (see, e.g., 
Ref. 2) and in many recent publications on MgB2.^'''^^"^^ 
An "operational" justification implied for this apparent 
contradiction usually goes like this: of course, the field 
distribution for a single vortex is described by Eq. (1) 
with a constant A. However, in the mixed state the aver- 
age order parameter is suppressed by overlapping vortex 
fields, and therefore, the "macroscopic" A (7^ to A cal- 
culated for H 0) should enter the free energy Eq. (2). 
This macroscopic parameter may depend on the average 
magnetic field B. 

The inconsistency of such an argument is exposed if 
one considers the clean limit at zero temperature. In 
this case the London A does not depend on the order 
parameter (in fact, it depends only on the total electron 
density), so that the mixed-state order parameter sup- 
pression cannot be reffered to as a general justification 
for employing A(i3). 
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D. On the superfluid density 

The quantity A^(0)/A^(r) is often taken as a measure 
of the superfluid density Ug. This assignment has unam- 
biguous justification only for isotropic superconductors 
(see, e.g., Ref. 30) and when A(T) is defined as the pene- 
tration depth of a smaU magnetic field (strictly speaking 
in the limit H ^ 0). One of the attractive features of the 
standard London theory is that one can extract 1/A^ (T) 
directly from the magnetization (3) by measuring the 
constant slope dM/d(\nB) for each temperature. This 
way of determining the superfluid density Us rests, there- 
fore, upon whether or not the standard London model 
for M is valid. As we have seen in a number of examples 
above, this is quite often not the case. 

Perhaps the best example of the futility of extracting 
rig from magnetization data is provided by the data for 
VsSi. As is seen in Fig. 9, e.g., for T = 14K, the slope 
dM/d{\nB) decreases with the field increasing up to ~ 
1 T, but it grows with further field increase toward Hc2- 
If one takes literally the proportionality between and 
the slope dM/d{\nB), one should conclude that the field 
suppresses the "superfluid density" as long as B is under 
« 1 T but for i3 > 1 T, Ug is enhanced by increasing field. 

Hence, extracting any quantitative information about 
the superfluid density or the penetration depth from the 
magnetization data with the help of the standard London 
model should be taken "with a grain of salt" at best even 
for such "simple" materials as VsSi, not to speak about 
MgB2 in which the different field behavior of the two gaps 
further complicates the matter. ^^'^^'^^ 



To summarize, the field dependence of the core 
size ^(B) has been incorporated in the London model 
to describe intermediate field reversible magnetization 
M{B,T) for materials with large k. The dependence 
^{B) is directly related to deviations in M{lnB) from lin- 
ear behavior prescribed by the standard London model. 
A method to extract ^{B) from the magnetization data 
is proposed. For most materials examined, £,{B) so ob- 
tained decreases with increasing field; the dependence 
becomes weaker with increasing temperature or with 
strengthening the non-magnetic scattering, in qualitative 
agreement with theoretical predictions and with existing 
/LtSR, SANS, and flux- flow resistivity data. The method, 
however, fails when applied to MgB2 and - surprisingly - 
to VsSi, the subject for a separate discussion. 
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